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Abstract We apply techniques of subriemannian geometry on Lie groups to laser-induced popula- 
tion transfer in a three-level quantum system. The aim is to induce transitions by two laser pulses, 
of arbitrary shape and frequency, minimizing the pulse energy. We prove that the Hamiltonian 
system given by the Pontryagin Maximum Principle is completely integrable, since this problem 
can be stated as a "k p problem" on a simple Lie group. Optimal trajectories and controls are 
exhausted. The main result is that optimal controls correspond to lasers that are "in resonance" . 
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1 Introduction 

1.1 Physical context 

In the recent past years, people started to approach the control of the Schrodinger equation, using techniques of 
geometric control theory (see for instance [|[ ^ |l9[ ^). In this paper we apply techniques of subriemannian 
geometry on Lie groups to the population transfer problem in a three-level quantum system driven by two 
external fields (in the rotating wave approximation) of arbitrary shape and frequency. The aim is to induce 
complete population transfer by minimizing the pulse energy. The dynamics is governed by the time dependent 
Schrodinger equation (in a system of units such that fi — V): 

^'-^-mt), (1) 

where V(-) : M ^ and: 

H = \ ni E2 Vl2 I . (2) 
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Here (*) indicates the complex conjugation involution. The controls ^2(-)i that we assume to be different 

from zero only in a fixed interval [0, T], are connected to the physical parameters by ^lj{t) = fj,jJ-j{t)/2, j = 1, 2, 
with J^j the external pulsed field and fj,j the couplings (intrinsic to the quantum system) that we have restricted 
to couple only levels j and j + 1 by pairs. 

Remark 1 This finite-dimensional problem can be thought as the reduction of an infinite-dimensional problem 
in the following way. We start with a Hamiltonian which is the sum of a "drift-term" Hq, plus a time dependent 
potential V{t) (the control term, i.e. the lasers). The drift term is assumed to be diagonal, with eigenvalues 
(energy levels) ... > E3 > E2 > Ei. Then in this spectral resolution of Hq, we assume the control term V{t) to 
couple only the energy levels Ei,E2 and E2,E^. The projected problem in the eigenspaces corresponding to 
Ei,E2, E3, is completely decoupled and is described by the Hamiltonian (||) . 

The problem is the following: 

Problem. As.sume that for time t < the .state of the .system lies in the eigenspace corresponding to the 
ground eigenvalue Ei. We want to determine suitable controls i — 1,2, such that for time t > T, the 

system reaches the eigenspace corresponding to E3, requiring that these controls minimize the cost (energy in 
the following): 

j= r {\ni{t)\^ + \n2{t)\') dt. (3) 

Jo 

In 1^ this problem was studied assuming that the controls Clj are "in resonance" : 

n,it) = u,it) c., = E,+^ - E,, 

Mj(.):K^M, aje[-TT,n], j = l,2. (4) 

In the sequel we call this second problem (of minimizing the cost (^, which in this case reduces to (^{ui(t)^+ 
M2(i)^) dt), the "real-resonant" problem. The first problem (with arbitrary complex controls) will be called 
the "general-complex" problem. 

In 1^, the real-resonant problem was treated as follows: 

• first, using a time dependent change of coordinates that leaves invariant the source (eigenstate 1) and 
the target (eigenstate 3), it was possible to eliminate the drift term and hence to reduce the problem to 
a singular- Riemannian problem over the real sphere S^; 

• second, the Hamiltonian system obtained from the Pontryagin Maximum Principle (PMP in the follow- 
ing) was Liouville integrable and explicit expressions for geodesies were found (although not so simple 
expressions). 

In this paper, we address both problems in a more abstract setting: 

• in both cases, first, we eliminate the drift. For the general-complex case, we use a time dependent change 
of coordinates plus a gauge transformation. Again, this change of coordinates leaves invariant both the 
source and the target; 

• second, in both cases, we lift the problem into a right-invariant subriemannian control problem on a real 
simple Lie group G, {G = G^ = 50(3) for the real-resonant problem and G — G^ — SU{3) for the 
general-complex one). This subriemannian problem has very special features: 

There is an element 50 of C G^, and subgroups K^, K^, (K^ « 0(2), K"^ « U{2)) such that, denoting 
by k*, the Lie algebras of K^, respectively: 

a. The Lie algebras L,k of the pair {G,K) have associated Cartan decomposition L = k © p, with 
the usual commutation relations, 

b. The right-invariant distribution is determined by the p subspace of L, 

c. The right invariant subriemannian metric is determined by a scalar multiple of Kil\^, where Kil\^ is 
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the Killing form restricted to p, 

d. The source S and the target T of the optimal control problem are respectively S = Kgg — gQKgQ~^, 
the conjugation of K by go, and T = g^S; 

• subriemannian problems on semi-simple Lie groups in which the distribution is determined by the p 
subspace of a Cartan decomposition and the metric is proportional to Kil\^ (called in the following 
(k © p) -problems) have two important features: 

— abnormal extremals (if any) are never optimal since the so called Goh condition (which is a necessary 
condition for optimality of abnormal extremals, see is never satisfied (see Appendix C); 

— normal extremals can be computed with a very powerful technique developped by Jurdjevic in 
]T^ , pj[ . In particular he proved that the Hamiltonian system given by the PMP is completely 
integrable and he gave explicit expressions for geodesies. For sake of completeness, the Jurdjevic's 
technique is summarized in Appendix B; 

• as the PMP requires, we apply suitable transversality conditions, corresponding to the source and the 
target. These conditions allow to restrict the set of admissible extremals. 

• we prove that all the geodesies satisfying the transversality conditions have the same length. Optimality 
follows. 



Remark 2 In the paper [|17| , the connection between the geodesies for the k © p problem and the geodesies of 
the Riemannian symmetric space G/K is made. But in our case: 

(a) we are not on the Riemannian symmetric space G/K, but on some other homogeneous space G/K, K C K , 

k^K- 

(b) as the reader can check, it is possible to project on the symmetric space G/Kg^, but this problem is not 

invariant by the action of G on G/Kg^. In fact, we get a singular Riemannian problem over G/Kg^ (as 
we have shown in the paper Q). 

With this methods, for the real-resonant problem we get in a more natural way the result of and we give 
much simpler expressions for geodesies and optimal controls. In the general-complex problem, besides finding 
explicitly expression for optimal trajectories and controls we prove our main result: 

Theorem 1 (Main Result) For the three-level problem with complex controls, optimality implies resonance. 
More precisely, controls ^2i-) o.fe optimal if and only if they have the following form: 

r f7i(t) = cos(t/V3)e'[(^=-^i)*+^il, , 
\ r!2(t) -sin(VV3)e'[(^^-^=)*+^^l. ^ ' 

where ^Pi,f2 a'^e two arbitrary phases. Here the final time T is fixed in such a way subriemannian geodesies 
are parameterized by arclength, and it is given by T — -^tt. 

This fact was pointed out as an open question in jsj. In other words, optimal trajectories for the real-resonant 
problem are optimal also for the general- complex one. 

Remark 3 The fact that optimality implies resonance was also proved in for the two-level case, reducing 
the problem to an isoperimetric problem. More precisely the two-level problem, which is described by the 
Hamiltonian: 

W{t) E2 J ■ 

can be reduced to a 3 dimensional contact subriemannian problem, with a special feature: it has a symmetry, 
transverse to the distribution. It is a standard fact that such a subriemannian problem is an isoperimetric 
problem (in the sense of the calculus of variations) on the quotient by the symmetry. In this case, the fact 
that optimality implies resonance is nothing but the (trivial) solution of the classical isoperimetric problem (or 
Dido problem) on the Riemannian sphere. We refer to 1^ for details. 
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Remark 4 Besides physical applications, the three-level problem is particularly interesting since it has a non- 
trivial geometric structure, but it is completely computable. On the other side the two-level problem, described 
in Remark ^ (although it is a k © p problem, with a nice geometric description) has a quite trivial solution. 

Systems with more than 3 levels appear to be more difficult to treat. The drift can be eliminated in the 
same way, but one gets a problem in which the control distribution is determined by a strict subspace of p 
only. In this case the Jurdjevic's formula fails to be true, and the integrability of the PMP is an open problem. 
See Remarks |l^ and ^ (the last in Appendix B). In this case the Hamiltonian reads: 



H 
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Remark 5 For questions related to invariance under time reparameterization of our optimal trajectories, we 
refer to 



1.2 A crucial fact 

At this level, we want to point out an obvious, but important symmetry property. It is important for compu- 
tations in the following, but also for practical applicability of the result. 

The optimal control problem we consider (in C'^ or in the complex sphere in C^) is invariant under mul- 
tiplication by a constant phase e''^. Hence, the problem of minimizing the energy from a fixed point in the 
ground eigenstate, to the third eigenstate, is the same as minimizing the energy to move from the full ground 
eigenstate, to the same target: actually, if {'il'{t),fl{t)) is an optimal pair trajectory-control, from -00 in the 
ground eigenstate to the third one, (e"^7/'(t), 57(t)) is an optimal trajectory from e^'^V'o- Moreover, as a conse- 
quence, whatever the initial condition iJjq in the ground eigenspace, the optimal control is the same . This is 
very important for applicability. Although it is a time dependent problem, these considerations hold also for 
the real-resonant problem. 



1.3 The problem downstairs and upstairs 

The problem of inducing a transition from the first to the third eigenstate, can be formulated, as usual, at the 
level of the wave function tp{t) but also at the level of the time evolution operator (the resolvent), denoted here 
byg(i): 

m^gimO), 9it)&Ui3), 9(0) ^td. (7) 

In the following we will call the optimal control problem for tpit) and for g{t) respectively the "problem down- 
stairs" and the "problem upstairs". 

The state- vector ip{t), solution of the time-dependent Schrodinger equation iip — H^, where H is given by 
formula (||) can be expanded in the canonical basis of C'^, formed by elements ipi — (1, 0, 0), ip2 — (0, 1, 0), ips = 
(0,0,1): V(t) = ci{t)ipi + C2it)ip2 + C3{t)ip3, with |ci(t)|2 + |c2(t)P + \c3{t)\^ = 1. For < < and t > T, |c.(i)|2 
is the probability of measuring energy Ei. Notice that, since rij(i) = 0, for alH < 0, t > T, j = 1, 2, we have: 

^|ci(i)p==0 for i < and t > T. 



At the level of the wave function, we formulate the problem in following way. Assuming |ci(t)| = 1 for <: < 



we want to determine suitable control functions j — 1, 2, such that |c3(i)p = 1 for time t > T, requiring 



that they minimize the cost (|3D . Thus we have a control problem on the real sphere C C'^ with initial point 

-rd 



be long ing to the circle defined by |cip = 1 and target defined by |c3p = 1. Equivalently, as we said 



in L2, the initial point V'(O) can be considered as free in S^. In the following the labels ( ) and ("), indicate 
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respectively downstairs and upstairs. Sources and targets upstairs (that will be called 5^ and 7^") will be 
computed in Section |4.3| after elimination of the drift. Why we use the subscript (c), will be made clear in 
Remark |^. 

1.4 Contents of the Paper 

The paper is organized as follows. In Section ||, we define precisely the k p problem and in Section ^, we 
discuss the elimination of the drift term. In Section ^ we formulate our problems in the k © p form and 
we define sources and the targets. In Section ^ we compute optimal trajectories reaching the final target and 
satisfying transversality conditions. The proof of Theorem |l| follows. In Appendix A we point out an interesting 
consequence of the Cartan decomposition, while in Appendix B and C respectively, we explain the Jurdjevic's 
formalism and why abnormal extremals are not optimal in a semi-simple k p problem. 



2 The k e p Problem 

For sake of simplicity in the exposition, all over the paper, when we deal about Lie groups and Lie algebras, 
we always consider that they are groups and algebras of matrices. 

Let L be a semi-simple Lie algebra and let us denote the Killing form by Kil{., .), Kil{X, Y) ~ Tr{adx°adY)- 
In the following we recall what we mean by a Cartan decomposition of L. 

Definition 1 A Cartan decomposition of a semi- simple Lie algebra L is any decomposition of the form: 

L = k®p, iw/iere [k, k] C k, [p,p]Ck, [k, p] C p. (8) 
Remark 6 Since L is semi-simple then relations (||) implies [p, p] — k, [k, p] = p (see Appendix A for the 
proof). This fact will be crucial for the elimination of abnormal extremals. 

Definition 2 The right- invariant k p control problem on a compact semi-simple Lie group G is the subrie- 
mannian problem with right-invariant distribution induced by p and cost: 

< 9 9^^, 9 9^^ > dt, 

Jo 

where < , >:= —a Kil\^{ , ), a > 0. 

The constant a is clearly not relevant. It will be used just to obtain good normalizations. 

In the following we will be interested in k® p problems on so(3) and su(3), for which we have respectively 
Tr{adx ° ady) — Tr{X Y) and Tr{adx o ady) = 5 Tr{X Y). Then, in order to get for both Lie algebras the 
useful relation: 

<X,Y>=-h'r{XY), x,yep, (9) 

we must set a = 1/2 for so(3) and a = 1/10 for su(3). 

Let {Xj} be an orthonormal (right-invariant) frame for the subspace p of L, with respect to the metric 
defined in Definition ||. Then the k p problem reads: 

g ~ "i^i 9^ min / dt, uj G R. (10) 

From relations (||), one gets that the so called "Goh condition" is never satisfied (see Appendix C). As a 
consequence: 

Proposition 1 In the problem defined in Definition^, every strict abnormal extremal (if any) is not optimal. 

Remark 7 The k © p problem can be also stated in the case where L is of noncompact type. In this case, to 
have a positive definite metric, we have to require that k is a maximal compact subalgebra of L, and we must 
define < , >:= -\-a Kill ( , ), a > 0. 
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3 Elimination of the Drift Term 

In this section, we show how to ehminate the drift term from a n-level system of the form with general- 
complex and real-resonant controls. For G C, let us denote by Mj(n) and Nj(n) the n x n matrices: 

where S is the Kronecker symbol: Sij — 1 if i — j, Sij — if i ^ j. Let A = diag{Ei, i?„), ujj = -Ej+i — Ej, 
j = — 1. We will consider successively the general complex problem (in dimension n): 

ri-l 

iip = Hip, H ^ A + Y^ where VI j e C, 



(this is nothing but another notation for the matrix (|^)), and the real-resonant problem: 

n-l 

H = A + ^Mj (e*("^- ) ) , Uj , a, e K. 

In both problems, 'ip lies in the complex sphere in C", and we want to connect the source 5{? = {(e*''', 0, 0)} 
to the target — {(0, 0, e*"^)}, by minimizing: 

n— 1 n — 1 

J = I 'y^\^j^dt, which in the real-resonant case is J = / 'S^u^^dt. 

Jo ^.^1 Jo ^.^1 

In both cases, we first make the change of variable ip — e~'^*A (interaction representation), to get (here 

Let us stress that the source and the target are preserved by this first change of coordinates. 
3.1 The general-complex case 

In that case, we make the time-dependent gauge transformation (i.e. cost preserving change of controls): 
Hence our problem become (after the change of notation A ip, flj —^Uj): 

fa. min/„^E;CiKNi, ^(0)e5^, x(r)erc^ 

Notice that the matrices Nj{l), Nj{i) generate su(3) as a Lie algebra. For n = 3, the Schrodinger equation 
(|l2|b) writes, in matrix form: 

/ ui(t) \ 

tP = Hc'iP where He \ -ul{t) U2{t) . (13) 

V -ul{t) I 

The cost and the relation between controls before and after elimination of the drift are: 

j = {\um^ + \u2{i)\^) dt, (14) 

f7i(<)=zii(t)e«[(^-^^)*+'^/2l, 

r!2(<) -U2(t)e'l(^^-^^)*+"/'l, ^ ^ 
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3.2 The real-resonant case 

In this case, since ilj = Uje**^'^^*"'""^\ we have: 

n-1 

iA = Y^ Mj (e*"^Uj)A, Uj G M. (16) 
i=i 

We make another diagonal, Hnear change of coordinates: 

A = e'^4>, L = diag{Xi, ...A„), for Ai, A„ € M. 

This gives: 

n-1 

Choosing the A^'s for e^^°'^'^^^+^~^^^ — i, we get: 

n-1 

The source and the target are also preserved by this change of coordinates. Notice that the matrices Nj{l),j = 
1 . . . n — 1 in (^), generate so{n) as a Lie algebra in its presentation by skew-symmetric matrices. This means 
that the orbit of the system ( p7| ) through the points (±1,0,... ,0) is the real sphere S*""^. In other words, the 
action on C" of the subgroup SO{n) C SU{n), restricts to the reals. Hence (by multiplication on the right by 
e**^), the orbit through the points ...,0) is the set S'"~^e*'''. Therefore, (after the change of notation 

(f) ^ ip) the real-resonant problem reduces to the problem over S"""^ (see sections for more details): 

min/o^E'Ci "'rfi, x(0)e{(±l,0,...,0)}, x(r) e {(0, 0, ±1)}, 

For n — Z we get for the Schrodinger equation in matrix form: 

/ ui(<) \ 
i) = HKip where #r := -Ui(t) U2{t) , (18) 

V -U2{t) ) 

The cost is given again by formula (^4[) and the relation between controls before and after elimination of the 
drift is: 

VLjit) = Uj{t) e*('^j*+"^-), ojj = Ej+i - Ej, 

Uj{.):R^R, aje[-TT,n], j = l,2. (19) 

Remark 8 In the following, besides to the labels {'^) and (") that indicate respectively downstairs and upstairs, 
we will use the labels (c) and (r) to indicate respectively the general-complex problem and the real-resonant 
one. When these labels are dropped in a formula, we mean that it is valid for both the real-resonant and the 
general-complex problem. With this notation: 

5^ = {(e^^0,0)}, Tc''-{(0,0,e'^)}, 
5^ = {(±1,0,0)}, rK'' = {(0,0,±l)}. 
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4 The k0p Form of the Problem, ControllabiUty, Sources and Tar- 
gets 

4.1 Controllability 

First notice that, after elimination of the drift, the problem upstairs is in SU{3) and not in J7(3) (the center 
of C/(3) is eliminated). Moreover, as we already observed, in the real resonant case, starting from (1,0,0), we 
have the standard action of SO{3) on M'^. From the fact that the matrices i/u generate so(3) as a Lie algebra, 
it follows easily that the orbit through (1, 0, 0) is the real sphere of equation: 

cl + cl + cl = l, CjGM, j = 1,2,3. (20) 

It is equivalent to say that the real-resonant problem is not controllable on S^. This fact is explained in 
details in Moreover in |^ it is proved that on each of these spheres 5^ the control problem reduces to a 
singular-Riemannian problem. The "relevant locus", which is the union of all the orbits (translations of S'^) 
passing through the eigenstate number 1, has an interesting non trivial geometric description. It is the only 
non orientable sphere-bundle over S^. 

After elimination of the drift, the real-resonant problem upstairs is on 50(3), and the general-complex 
problem upstairs is in SU{3). We have: 

Proposition 2 After elimination of the drift, the real resonant problem upstairs is completely controllable on 
Gr := SO{3) C U{3), and the general- complex problem is completely controllable on Gc '■— SU{3) C U{3). 
As a consequence, by transitivity of the action on the real and complex spheres, the corresponding problems 
downstairs are controllable on S"^ and . 

Proof. This is a consequence of the general theorems of controllability of left-invariant control systems on 



compact groups: the Lie-rank condition is necessary and sufficient for controllability (see |18 ) 



Remark 9 In fact, before the elimination of the drift, the general complex problem upstairs is in C/(3), due to 
the nonzero trace of A. For the same reason as in Proposition]^, it is completely controllable. 

4.2 The problems in the k © p form 

For the lifted problem, from ?/> = H^, using (Q) one gets: 

9 = Hg. (21) 

where generates so(3) and He generates su(3) (as Lie algebras, for distinct values of the controls). For 
both Hamiltonians, the problem of minimizing the cost ( p^ is a k © p problem. 
In fact in the real-resonant case equation (pi|) can be written as: 



g = [uiXi + U2X2)g, 



where: 





1 


















• 




• 





■) 













-1 





Setting: 



= -1 , ^2 = 1 . (22) 



X3 =( |, p : =span({Xi,X2}), and k = span({X3}), 
-10 
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one gets relations (^) (with L = so(3)). Moreover the distribution is right-invariant and the frame ( p2| ) is 
orthonormal for the metric (||). 
In the case of sm(3) we have: 

g = {uiXi + U2X2 + U'iYi + UiY2)g^ uj £ R, 

where Xi and X2 are given by formula and Yi, Y2 by: 





i 






( 







: 





• 










: 











^ 


i 





One can easily check that relations (||) are verified with p :=span({Xi, X2, Yi, I2}) and k :=span({Zi, Z2, Z3, Z4}), 
where: 



1 
Zi=\ I , Z2 
-10 



i 

I , Z3 

1 





are the remaining 4 generators of sm(3). Again the distribution is right-invariant and the frame {^1, X2, Yi, I2} 
orthonormal for the metric (j^). 

Remark 10 • (general complex case) k is a subalgebra of su{3) containing a maximal Abclian subalgebra of 
su(3), which is generated by Z^ and Z4. From (see also [^) it follows that k must be a Borel subalgebra. 
In this case, it is u(l) x su{2) — u{2). 

Let us denote by S{U{1) x C/(2)) (resp. S{Z2 x 0(2)) , the groups of matrices of the form: 





( e 


^ 




B = 





J 


, where det(i3) — 1, 



with e e C/(l), (resp. e e Z2 ^ {-1,1}), and U, £ U{2) (resp. t/, e 0(2)), S'(C/(1) x [/(2)) w C/(2), 
^(^2 X 0(2)) « 0(2). 




Set: 

ffo := I 1 I , (23) 
V 1 / 

go e S'0(3) C SU{3), go^ = 50^^- In the real-resonant case, we will denote by the subgroup of — 50(3), 
conjugate by go~^ to 5(^2 x 0(2)) « 0(2), with Lie algebra k^ = {Xs}^^. 

In the general-complex case, we denote by K'^ the subgroup of = SU{3) conjugate by 50^^ to S{U{1) x 
C/(2)) « C/(2) with Lie algebra k^ = {Zi, Z2, Z3, Z^}^^. 

• In the general-complex case, the corresponding Riemannian symmetric space G/K is SU{3)/ S{U{2) x 
U{1)) ~ CP^ that has rank 1 (see Helgason pp. 518), and the Cartan subalgebra of p has dimension 
1 only, as one can check easily. 

• In the real- resonant case, it is RP'^, the 2-dimensional projective space. 

Remark 11 Notice that in the real- resonant problem the distribution is a contact distribution. In this case it 
is a standard fact (see for instance ^ ^ ^ ^) that there are no abnormal extremals. On the other hand, 
the distribution of the general complex problem is not a contact distribution and abnormal extremals do exist. 
Anyway they are not optimal, due to the fact that the ''Goh condition" (see Appendix C) fails to hold. 
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Remark 12 (n-level case) For n > 4 the Hamiltonian given in formula (12) generates so{n) or su{n) respec- 
tively with real or complex controls . Anyway it never gives rise to a k©p problem since the distribution is only 
a strict subspace of p. As explained in Remark ^ this fact causes the failure of the proof of the integrability 
of the PMP. No results in this paper about optimal trajectories can be easily generalized to 4 or more levels. 
For 4 or more levels the integrability of the PMP, and the fact that optimality implies resonance, are open 
questions. 



4.3 Sources and Targets 

In this section, we describe sources and targets for the real-resonant and the general-complex problems upstairs. 
As we already said, from section 1.2, sources and targets downstairs, for the general-complex problem and 



the real-resonant one, are respectively the circles S^, T^, and the sets 5^ — {(±1,0,0)}, 7^ ~ {(0,0, ±1)}. 

In both cases, we decide that the canonical projection tt : G — ^ G/K maps K to the point (1,0,0) in the 
complex and real sphere in C'^ (resp. M'^). Here Kc — SU{2) (general complex case) and Kr — SO{2) (real 
resonant case). Let 5" = n-^S"^). Then = S{U{1) x C/(2)), 5| = 5(^2 x 0(2)). 

The element go G 50(3) defined in the previous section maps to T^. Then T" :— 7r^^(T'^) — goS"" in both 
cases as it is easy to check. 

It is clear that a trajectory g{t) of the system upstairs: 

g = Hg, te[0,T] 

such that g{0) e 5" , g{T) e T", maps into a trajectory g{t)K of the system on the sphere, 

gK = Hgk g{0)K € S'^, g(T)K G T'^ 

with the same control, hence the same cost. Conversely, if x{t), t e [0, T], is a trajectory on the sphere, of the 
system i = Hx, that maps the point (£,0,0), |e| = 1 to the point (0,0, e'), \e'\ = 1. Then, the corresponding 
fundamental matrix solution g{t) satisfies, for all s e 5", g{0)s ~ s E S"^ : 

ff(r).s.(£,0,0) = (0,0,£')=5o.(£',0,0) 

Therefore, g^'^g{T).s.{e, 0,0) = (e',0,0) and g^^.g{T).s G 5", g{T) e go^" = T"". This shows that, for all 
s E 5", the lifted solution upstairs satisfies g{0).s G 5", g{T).s e T^. 
Hence we get the following table: 



PROBLEM 



SOURCE 



TARGET 



real-resonant 



Z2 








0(2) 



e 50(3) ) = S{Z2 X 0(2)) 



:= go 5" = goS{Z2 x 0(2)) 



general-complex 



5" 



C/(l) 








U(2) 



e SC/(3) \ = S{U(1) X t/(2)) 



r« := goS^ = goS{U{l) X U{2)) 



5 Expression for Geodesies, Controls and Transversality Conditions 
5.1 Preliminaries 

As explained in Appendices B and C, candidates optimal trajectories for the k p problem are only "normal 
geodesies" and are given by the following formula: 



m = e 



— p-'^kt {Ak+Ap)t 



5(0), 



(24) 
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where 5(0) is the starting point belonging to the source. Set M(0) = Mp{0) + Mk{Q) = A = Ak + Ap E L for 
the initial value of the covector P(0) — d*^ Pg{a) — {M{0), •), where (•, •) is given in (^), and with the notations 
of Appendix B. Here, Ap = Mp{0), Ak — Mk{0) are the orthogonal projections of M{0) = A on p and k 
respectively. Geodesies are parametrized by arclength iff: 

{Ap,Ap) = l (25) 

5.2 Transversality conditions 

In both cases, our source is Kg^ = cjoKgo^^ where go has been defined in formula ( p3| ) and K'^^ = S{U{1) x U{2)), 
Kf^ = S{Z{2) X 0(2)). Let us notice the two following facts: 

Facts: 



1. Transversality conditions at the source may be required at the identity only. 

2. Transversality conditions at the source imply transversality conditions at the target. 

Proof : Point 1 comes from right-invariance and the fact that the source is a subgroup. Point 2 comes from 
the following lemma |l|. H 

Let t E [0,T] be a normal extremal corresponding to the covector Pg[t), such that g{0) = Id. 

Lemma 1 pidiAdgak.) = implies pgi^i^idL g(^t)Adg„'k) ~0,\/t E [0,T]. 

Proof : (with the notations of Appendix B). Set / = Pg(^t){dLg(^t'jAdgg'k) = dR*^^^_iP{t){dLg(t)Adgg'k). Then 

/ = P{t){Adg^,^g„k) = {M{t),Adg(^t)g„k). But: 

M{t) = Mpit)+Mk{t) = Mp{t) + Mk{0), 

M{t) = e-^'^^(°)*Afp(0)e*^'=(")* + A4(0) = e-^^'=(°)*(Mp(0) + Mfc(0))e*^^(")* = Ad^-M,mtM{0), 
and g{t) ^ e-^^^(")*e*^(o)*. Then: 

/ — {Ad^-Mi^mtM{0), Ad^-M^mtAd^Mwt Adg^k), 
and since the killing form is Ada invariant, / — (Af(0), Ad^M(o)tAdg„]i.). Hence, for the same reason: 
/ - {Ad^-M(oHM{0),Adg„k) - {M{0),Adg„k) = P{0){Adg„k) - pjdiAdg^k) = 0. 

■ 

Similarly to A^i, N2, (cfr. formula (|ll|)) let us define by: 

/ ase^^^ \ 

NiAase'"') = 0. (26) 

\ -age-'^a 0/ 

Let us set, in the real-resonant case Ap ~ aiiVi(l) + a2N2{l) and Ak — a3iVi,3(l). In the general complex case, 
set Ap = iVi(aie^^i) + 7V2(a2e'^"), Ak = 04^3 + 05^4 + A^i,3(a3e'^3)^ Here, a, G R and E [~n,n]. 

5.3 The real-resonant case (G=SO(3)) 

Proposition 3 For the real-resonant problem, the transversality conditions Kil{A, TidS^) — implies 02 = 0. 
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Proof. We have: 



AdgX 



rp (PU 





























/5e 



Then equation Kil{A, Ti^S^) = is satisfied for every /3 e R if and only if 02 = 0. 
From Proposition ^ and condition (^5|), one gets the covectors to be used in formula 



±1 03 
Tl 
-as 



(27) 



Proposition 4 The geodesies ({24i), for which g{0) = Id, with A given by formula \2'\), reach the target 
for the smallest time (arclength) \t\, if and only if — ±l/-\/3. Moreover, the 4 geodesies (corresponding A^ 
and to the signs ± in a^) have the same length and reach the target at time: 

T ^ 

1 = TT. 

2 

Proof. Computing g{t) = Q-Akt^{Ak+Ap)t ^ with A given by formula (|27|), and recalling that: 

m = mm = |^ o j , 

one gets for the square of the third component of the wave function: 

2 

(cos(ta3) sin(t vT+ap^) 03 vTTos^ — cos(t vT+a?) sin(ta3) (l + as^)) 



(C3(i)) = 



By lemma (2) in Appendix D, we get the result. 



[l + as^f 



(28) 



Explicit expressions for the wave function and for optimal controls 

Let us fix for instance the sign — in ( |27| ) and 03 = +l/\/3. The expressions of the three components of the 
wave function are: 



ci(<) =cos(^)^ 
c,{t) = ^ sin(^) 
csit) = -sm{^f 

Let us stress that this curve is not a circle on S'^. 
Controls can be obtained with the following expressions: 



(29) 



(30) 



We get: 



ui{t) = -cos(-^) 
U2(<) =sin(^) 



(31) 



The situation is depicted on Figure 1. 
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Figure 1 

The result given by formulas ( |29| ) and ( pT| ) is the same result as the one of but it has a simpler form. Using 
expression (|9h (resonance hypothesis), we get for the external fields: 



ni{t) = -cos(V\/3)e*("i*+"i), 



(32) 



Notice that the phases ai, a-i are arbitrary. 

5.4 The general-complex case: G=SU(3) 

Proposition 5 For the general- complex problem, the transversality conditions Kil{A,TidS'^) — implies ai 
04 = 05 = 0. 

Proof. We have: 

Then equation Kil{A, TidS^) = 0, is satisfied for every ai, a2, Pi, (32 G K if and only if 02 = 04 = 05 = 0. 
The covector to be used in formula (p3) is then: 



ioii 











i{a2 — 


ai) Pi +1(32 





+ 


iP2 —iOL2 



e 

-e-*''i 
-036-*^^ 



(33) 



Proposition 6 The geodesies with A given by formula ( pj^ (for which g{0) = Id), reach the target 

for the smallest time (arclength) \t\, if and only of = ±l/\/3. Moreover all the geodesies of the two parameter 
family corresponding to 9i,6y, G [— tt, tt], have the same length: 

1 — TT. 

2 



Proof. The explicit expression for |c3p is given by the right-hand side of formula (28). The conclusion follows 
as in the proof of Proposition (^. 
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Explicit expressions for the wave function and for optimal controls 

The expressions of the three components of the wave function and of optimal controls are: 

ci(i) =cos(^)^ 

C2(t) = -f sin(J|)e-«^ (34) 
C3(t) = -sin(^)V^«3, 

r = cos(VV3)e*«i , . 

\ ^i2(^) = -sin(^/^/3)e^(^3-eO, ^'^^) 



Remark 13 Again, notice that none of these trajectories is a circle on the corresponding (translation of) sphere 
S^. Notice that all the geodesies of the family described by Proposition (^ have the same length as the 4 
geodesies described by Proposition ^. This proves that the use of the complex Hamiltonian (13) (instead of the 
real one (|l^)) does not allow to reduce the cost (p^, and this proves Theorem |l| in the introduction. 

Formulas (|9|), (|l]) can be obtained from formulas (|34|), ( |35| ) setting 6i = tt, 6:^ = 0. 

For the general-complex problem, using expressions (^5|) in ( p^ ) one gets for the external fields: 

f^i(i) = cos(t/V3)e'("^*+''i+f), 

n^it) = - sin(t/V3)e*("=*+«3-»i + f », ^ ' 

where uji = E2 — Ei, 0^2 = -E3 — -B2 ■ 

Formula ( ^6| ) coincides with formula (|^) setting ipi :— 61 + tt/2, ipi :— 9^ — 9i — tt/2. 

We recall that in previous papers for the real-resonant problem, (see Q for more details), optimality was set 
as an assumption while here, it is obtained as a consequence of the PMP. 

Moreover we get that the optimal cost and the probabilities |cj(t)p, j — 1,2,3, are independent from 
6'i,6l3 e [-7r,7r]. 



Appendix A: An interesting consequence of the Cartan Decomposition 

Theorem 2 Let 1, be a semi-simple Lie algebra and L = k©p a Cartan decomposition (i.e. it holds [k, k] C k, 
[p, p] C k, [k, p] C p). Then [k, p] = p, [p, p] = k. 

Proof. We can restrict ourself to the case in which L is simple, since a semi-simple Lie algebra is a direct sum 
of simple ideals. 

Claim 1: p + [p, p] is an ideal in L 
Proof of Claim 1: 

• [Pi P + [Pi P]] — [Pi P] + [Pi [P) P]]- Using the Cartan relations it follows that the second term is contained 
in p. Hence [p, p + [p, p]] C p + [p, p] . 

• [k, p + [p, p]] = [k, p] + [k, [p, p]]. Now, using the Cartan relations we have: [k, p] C p while [k, [p, p]] = 
— [p, [k,p]] — [p, [p,k]] C [p,p], where we have used the Jacobi identity. Therefore [k, p-|-[p, p]] Cp-|-[p,p]. 

Claim 2: k + [k, p] is an ideal in L 

Proof of Claim 2: 

• [p, k + [k, p]] = [p, k] + [p, [k, p]] C [k, p] + k. 
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. [k, k + [k, p]] = [k, k] + [k, [k, p]] C k + [k, p] . 

From the fact that L is smiple (the only ideals are and l) it follows p + [p,p] = L [PjP] — k, 
k + [k, p] = L ^ [k, p] = p. ■ 

Notice that in general is false that [k, k] = k. In particular it is false for the Cartan decomposition of SU (3) 
we used in this paper (k has a center). Theorem |^ is crucial to prove that, in our case, abnormal extremals are 
never optimal (see Appendix C). 



Appendix B: The K + P Problem 



Here following and , we recall how to write the PMP for a right-invariant control system on a Lie group 
G. We do it for a group of matrices only. 

Let L and L* be respectively the tangent and cotangent planes at the identity of G. Consider the right- 
invariant system: 

9 = X{u)g, (37) 

where X{u) E L and u belongs to the set of values of controls U. To have a right-invariant optimal control 
problem, we must also assume a right-invariant cost (that is a cost that does not depend on g, the coordinate 
on the group): 

/ f{u{t)) dt, where / : [/ ^ M is a smooth function. (38) 
Jo 

Moreover we assume that the initial and final points belong to two given smooth manifolds: 

9{0)eAhn, .g(T)eM/,„. (39) 

For each A S M and u £ U define the Hamiltonian : 

Hu{.) : L* ^ M, Hu{P) := P{X{u)) + Xf{u). 

We have dHu{P) G (l*)*, (P G L*) and since L has finite dimension, we can identify (l*)* ~ L and consider 
dHu{P) S L. Let B G L and denote with ad*g{.) the operator from L* L* defined by: 

[adUP)KA) :- P{adB{A)) = P([B, A]), V A e L. 

The PMP (see [^) is a necessary condition for optimality. For right-invariant control systems it reads (see 

Theorem 3 (PMP for right-invariant control systems) Consider a couple {u{.),g{.)) : [0,T] U x G 
subjected to the dynamics ([g^ and to the constraints (^^. // it minimizes the cost (f?^, then there exists a 
constant A < and a never vanishing absolutely continuous function P(.) : t G [0, T] i— > P{t) G L* such that: 

'^ = dH^(,){P{t))g{t) (40) 
^--«rfS«„.,(P(*))™)- (41) 



Moreover: 



Hu(t){P{t))='H{P{t)) where niP{t)) := m£ixR,{P{t)). (42) 
P(0)(Tg(o)Mj„.g"^(0)) = 0, P{t){Tg^T)Mf^n.g^\T)) = (transversality conditions) (43) 
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Here g{t) £ G and P{t) £ L* is the covector translated back to the identity. P{t) is related to the usual covector 
Pait) e T*i^t)G by: P{t) = dR^^^^pg^t) or pg^t) = dR^^_,Pg(t) 

The couples {u{.),g{.)) satisfying conditions (40), ( |4l| ) and ( ^2| ) with A = are called abnormal extremals. 
Couples corresponding to A 7^ (in this case we can normalize A = —1/2) are called normal 

extremals. 

In the following we show that for normal extremals this Hamiltonian system becomes completely integrable 
if the problem is k © p. Abnormal extremals may exist in general, but they are never optimal as explained in 
Appendix C. 

Consider the right-invariant k0p problem defined in Definitions ^ ||. Here L — k®p can also be noncompact 
(in that case k must be the maximal compact subalgebra of L, cfr. Remark 

The Killing form defines a non-degenerate pseudo scalar product on L. This permits to identify L with L* 

by: 

P e L* < — > M e L P(C) = Kil{M, C), y C £ L. (44) 
Let us "translate" equation {nW for M G L. For each C G L we have: 



KU{^,C) ^ '^{C)^-[ad*,^^p^,y^{Pm{C)^-P{t){[dH{P{t))M) 
= -Kil{M{t),[dH{P(t)),C]) ^ +Kil{[dH{P{t)),M],C), 

where we have used the invariance of the Killing form under Lie Brackets: Kil{[A, B],C) — Kil{A, [B,C]), 
that can be easily checked. The equation for M{t) is then in the famous Lax-Poincare form: 

= [dH{P{t)), M{t)], M{t), dH{P{t)) £ L. (45) 

Let {Xj} be an orthonormal (right-invariant) frame for the p part of L, with respect to the metric defined in 

Definition ^. We have then X{u) = UjXj and f{u) — u\ + ... + u^^ (here u = (ui, U2, Unp) and Up is the 

dimension of the p subspace) . Moreover decompose M = Mp + Mk , where Mp £ p and Mk £ k. 

Proposition 7 We have dH{P(t)) = X{u{t)) = Mp{t). 

Proof. The first equality can be obtained comparing equation g{t) = X(u(t))g{t) with equation ( ^0| ) and using 
the fact that the differential of the right translation is a linear isomorphism. Let us consider the second one. 
From the maximum condition (^) with A = —1/2 and f{u) = u\ + ... + uf^^, we get: 

u^{t) ^ Pg^t){t)iX,g{t)) = P{t){X,) ^ Kil{M{t),Xi) = Kil{Mp{t),Xi) because X, = X,{e) £ p. 

In the last equality we used the fact that Mp and are orthogonal subspaces for the Killing form. Therefore 
Mp{t) = E, Ktl{Mpit),X,)X, = u,{t)X, = Xiuit)). 



With Proposition the equation for M become: 

^ + ^ = [Mp{t),Mp{t) + Mk{t)] = [Mpit), Mk{t)]. (46) 

Using the Cartan commutations relations (||), we have [Mp{t), Mk{t)] C p and equation ( ^ splits into: 

dMk 



dt 



O^Mfe(t) = Affe(O), 



^ = Wlp{t),M,m. 

Hence all the k-components of the covector are constants of the motion. Integrating the equation for Mp, and 
setting Ap ~ Mp{0), Au := A4(0), we get: 

Mp{t) = e-'^'='^pe^'=* = e-'"^'^"' Ap. 
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From equation with dH{P{t)) = Mp{t), we have: 



dt 

and the solution is equation (p4|): 

g{t) = e-^^*e(^^+^p)*5(0) (47) 
Setting A -.^ Ap + A)^, the transversahty conditions ( ]4^ ) reads in L: 

/^^^(AT;(o)A^^n.5"'(0)) - 0, (48) 
i^iZ(M(i),Tg(t)M/„.g-i(t)) = 0. (49) 

Remark 14 (n-level case) In the case in which the distribution is only a strict subspace of p, as in the n-level 
case {n > 4), dH{P{t)) = X{u(t), but it is not equal to Mp{t) in general. Then equation (^) become: 

dMr, dMu r , , ^ s , . N , 

^ + ^ = [X{u{t)), Mp{t) + Mfe(i)], 

where now the right-hand side is not completely contained in p. This means that not all the k-components of 
the covector are constants of the motion, and the proof of the integrability of the Hamiltonian system given 
above fails. 



Appendix C: The Goh Condition 

Let us consider any subriemannian metric over a manifold M, defined by its orthonormal frame {X^, i = 
1, completely nonintegrable. Then a necessary condition for a strictly abnormal extremal (i.e. an abnor- 

mal extremal which is not normal at the same time) to be optimal is that it satisfies the Goh condition. 

Definition-Theorem (Goh-condition) Let {p{t),x{t)), t G [0, T] be the Hamiltonian lift of the abnormal 
extremal x{t) (see Theorem^. Then: 

{p{t){x,{x{t)) = G, p{t)^o, te[o,r], (50) 

and a necessary condition for optimality of x{.) is {p{t){Xi),p{t)(Xj)}{x{t)) = 0, t G [0, T], or: 

p{t){[X,,X,]{x{t)))=0, ie[0,T]. (51) 

This Theorem is a consequence of a (highly non trivial) generalized Maslov index theory developped in |4| . For 
our right-invariant problem, the relations (^0|), (|5|) give: 

p{0){X,{2d)) = 0, p{0){[X,,X,]){id)) = 0, 

they imply with Theorem ||of Appendix A, that p(0)(L) — 0. Then p{Q) = 0. This is a contradiction since p{t) 
has to be nonzero, for all t. 

Hence, strictly abnormal trajectories are never optimal in our k © p problem. 



Appendix D: A technical computation 
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Lemma 2 Set fa = cos(ta) sin(f\/l + a^) 



Moreover, \fa\ = liff 



\a\ 



2k ^ k 



cos(t\/l + a?) sin(ta), then \fa\ < 1. 



< 1 ; fc 7^ and t = 



In particular, the smallest \t\ is obtained for k = ±1, a = t = ^^'^ 

Proof. Set A = 



= tVl + a^, then 

= A cos(A0) sm{e) - cos{e) sin(A(9) 

= ( ( A cos( A6I) , sin( A6') ) , (sin {9) , - cos{9))) 



Both vi, V2 have norm < 1 and \fa\ < 1. Hence, for \fa\ = 1, we must have ||f 1 1| = ||f2|| = 1) ^^i 
that cos(A^) = and cos(6') = ±1. Hence 6 = kw, Xe = ^ + k'n, A = ^ + ^. Therefore, ' 
Conversely, choose k,k' meeting this condition, and 6 = kn. Then cos(^) = ±1, sin(A^) = ±1, fa{t) = ±1. 



ztV2- It follows 

= A < 1. 



J_ , k^ 
2k ^ k 



Now, \t\ = ^^'^^2 i and the smallest \t\ is obtained for fc = ±1 (if fc = 0, ^ = and fa{t) = 0). Moreover, 
^ + Y < 1 is possible only for {k,k') = (1,0) or (1,-1) or (-1,0) or (-1,-1). 
In all cases, |A| — 



2' 



andt = ±^. 
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